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Abstract — This research explores the role and representation 
of networlt structure for LTI Systems. We demonstrate that 
transfer functions contain no structural information without 
more assumptions being made about the system, assumptions 
that we believe are unreasonable when dealing with truly com- 
plex systems. We then introduce Dynamical Structure Functions 
as an alternative, graphical-model based representation of LTI 
systems that contain both dynamical and structural information 
of the system. We use Dynamical Structure to prove necessary 
and sufficient conditions for estimating structure from data, 
and demonstrate, for example, the danger of attempting to use 
steady-state information to estimate network structure. 

I. Introduction 

One of the fundamental issues for modeling, identifying, 
and controlling complex networked systems is inferring 
system structure from input-output data. Structure is often 
the key for understanding a variety of complex systems 
because it enables a decomposition of the complete system 
into an interconnection of subsystems. When analysis of the 
subsystems is comparatively simple, and the interconnec- 
tion structure is well understood, then the behavior of the 
complex system can be deduced from an understanding of 
its components. Moreover, exploiting structural information 
can tremendously reduce the conservatism of robust solutions 
designed to compensate for system uncertainty. This impact 
on complexity and uncertainty makes structural information 
extremely important in the analysis of complex networked 
systems. 

Examples of scientists working on identifying or ex- 
ploiting network structure arise in a variety of disciplines. 
Social scientists have developed a rich literature on the use 
of network models to describe interpersonal associations, 
perhaps one of the most famous works being Milgram's 
"small world" experiment in the 1960's in which letters 
passed from person to person were able to reach a particular 
target individual in only about six steps [15]. More recently, 
attention has focused on networks of business communi- 
ties [5], [13], [16], internet-enabled virtual communities 
[10], citation networks in scientific communities [18], [20], 
preference networks for product recommender systems [7], 
[19], distribution networks [1], and the detection and desta- 
bilization of terrorist networks [3]. Epidemiologists have 
developed models for the dynamics of both epidemic and 
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endemic diseases spreading through population networks [9], 
computer scientists have developed algorithms for searching 
over networks that are deployed in a number of popular 
applications [2], and biologists use microarray and other 
data sources to infer the regulation structure in genomic, 
proteomic, and metabolic networks [8], [14], [21], [22] . 

Discovering structure from data, however, can be difficult. 
Typical identification methods do not emphasize structure 
estimation, but rather focus on behavior generalization by 
selecting a model that accurately predicts system outputs 
for unobserved inputs. As long as the dynamic behavior of 
the system is accurately described, the question of structure 
is often avoided altogether For many applications, various 
model structures for the same input-output map are equally 
useful for forecasting and control. Nevertheless, sometimes 
it is important not only to describe the system dynamics 
accurately, but to do so with a model that correctly represents 
the structure of the original system. 

In contrast with these identification methods that em- 
phasize system dynamics over structure, inference methods 
have been developed that emphasize structure over dynamics. 
These methods employ graphical models to describe network 
structure. Nodes represent system states, understood to be 
random variables, and edges indicate conditional dependence 
between variables. Using Bayes rule, measurements can then 
be used to update prior distributions believed to characterize 
relationships throughout the network. A rich literature has 
grown in this area, and even issues of inferring causality 
from correlation have been addressed at some level [11], 
[12], [17]. 

Nevertheless, although these Bayesian Networks provide 
an efficient way to parameterize the joint probability distribu- 
tion characterizing the entire system, conditional probabili- 
ties do not capture system dynamics, and the most successful 
inference techniques only work on directed acyclic graphs 
[4]. For some applications, such as modeling the citation 
network for a particular body of research, assuming the 
network is acyclic is reasonable since papers generally only 
cite previously published work. There are many applications, 
however, such as modeling biological or social or economic 
networks, where such an assumption insisting on the absence 
of feedback dependencies between system states would be 



entirely unreasonable. Moreover, often an accurate represen- 
tation of system dynamics is as important as that of system 
structure. In these situations, new methods are needed. 

This paper introduces Dynamical Structure Functions as 
a structurally accurate representation of complex LTI sys- 
tems that do not ignore system dynamics. We begin in 
the next section by demonstrating that transfer functions 
contain no structural information without more assumptions 
being made about the system, assumptions that we believe 
are unreasonable when dealing with truly complex systems. 
We also highlight some common pitfalls when estimating 
structure from data. We then introduce in Section [ITJ the 
Dynamical Structure Function of an LTI system and discuss 
its properties. Section Hvl then uses Dynamical Structure to 
provide necessary and sufficient conditions for estimating 
structure from data, and an example is provided illustrating 
the danger of using only steady-state information to estimate 
structure. Section |V] then concludes with a discussion of 
future work. 

II. Background: Structure Estimation and 
Dynamic Systems 

Consider the network characterized by the linear system 



C = 



X — Ax + Bu 
y ~ Cx 



(1) 



where xeR",ue R"\ yeR'\ and C = [/ 0] . We are interested 
in inferring the causal dependencies between the p measured 
states, y, from limited data. Typically, m < n, p < n, and « 
itself is unknown. 

In this work we do not assume that the system Q is both 
controllable and observable from the particular inputs and 
outputs specified by u and y. In the complex systems context, 
such an assumption would be unreasonable to impose since 
the number of inputs and outputs is assumed to be very 
small compared to the total number of states. Thus, assuming 
controllabiUty and observability would be restricting our 
attention to networks with very special structure. As a result, 
we can not assume that Q is a minimal realization of the 
corresponding input-output transfer function, G, given by 



Gis)^C{sI-A 
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In this work we also do not assume knowledge of the sys- 
tem's order. Thus, the true system, Q, has a particular causal 
structure and complexity that we can only detect through 
our interaction with the system at u and y. Nevertheless, we 
do assume throughout this work that the transfer function, 
G, can be obtained from the available data, u and y, using 
standard identification methods. 

Notice that the transfer function does not directly reveal 
structural information of the system. For example, consider 
the system 
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Note that this system has a very definite ring structure, where 
xi — >X2 — > xj — > X4 — >xi. Nevertheless, the associated 
transfer function, G, is given by 
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where p{s) = s"^ + + + As + .9375, which reveals 
nothing about the ring structure of the system. Although the 
structure is easy to read from the actual state realization of 
the system, a transfer function identified from input-output 
data-even ;/ identified perfectly-does not directly yield any 
structural information about the system. 

Given this difficulty using the transfer function to obtain 
structural information, one may ask why not identify the 
state space realization directly. Nevertheless, it is difficult 
to identify a realization of the system without knowing the 
order of the system. In this work, we assume that structural 
information must be obtained from limited data, that is, 
with measurements that constitute only part of the complete 
state vector Moreover, we do not assume knowledge of 
the full system complexity, or true system order Later, we 
demonstrate how incorrectly assuming knowledge of the 
system order can lead to erroneous structural estimates. 

Thus, transfer functions are generally obtainable from 
input-output data, but they contain no structural information. 
At the other extreme, state space realizations contain all 
information about the system, but they are difficult to obtain 
from limited information. We are interested in something in 
between, a representation that may still be obtainable from 
input-output data, but that also contains information about 
both the dynamics and the structure of the system. 

Structure is typically represented by a graph. Nodes 
represent system variables, and edges represent interaction 
between variables. Directed edges capture notions of directed 
influence, often quantified by conditional probabilities. We 
will consider a directed edge to indicate a causal relationship 
between variables. Regardless of how the notion of directed 
influence is represented, however, the absence of an edge 
between variables indicates a kind of independence between 
those variables; zi — > zi instead of zi ^ zi means that zi 
does not depend directly on Z2- That is, any influence Z2 may 
have on zi may only occur indirectly through Z2's influence 
on other explicit variables (nodes) in the network, and their 
direct influence, in turn, on zi- In particular, it is critical 
to note that z\ ^ zi means there may not be some hidden 
variable, z,, that has not been represented in the graphical 
network model through which Z2 influences z\ - For structure 
to have meaning, even hidden, unmodeled variables should 
respect the graph defining the network and only operate 
within edges. This has important implications for dealing 
with uncertainty. 

In its simplest form, then, structure is simply a square 
binary matrix S with stj — 1 indicating the presence of an 



edge directed from zj to z,. For the system T given in Q 
we would define our explicitly modeled variables to be z = 
[zl; z2][y; u]. Simple structure, S, would then be a p + m by 
p + m binary matrix; for the example (|3} we would have: 
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Note that we consider that variables may automatically in- 
fluence themselves since they may be recursively generated, 
thus the diagonal of S is identity. The blocks Qt, Pt, Qf, and 
Pp correspond to the partition of z as inputs and outputs of 
the system T. The input structure, Pt describes how inputs, 
M, influence the measured variables, y. The output structure, 
Qt describes how the measured variables, y, influence each 
other. Under the interpretation that y corresponds to part of 
the system state vector, the output structure Qt may also 
be called the internal structure of the system T (with Pt 
then being called the control structure of T). The remaining 
blocks, Qf and Pf, describe the feedback environment of 
T . In general, when T is in feedback with an operator F, 
Pf is the input or control network of the feedback operator 
F, while Qf is F's internal or output structure. When no 
feedback operator is defined and the inputs u are truly 
considered free variables, then Pf is zero (since u does not 
depend on y), and Qf is identity (since m's only depend on 
themselves and do not influence each other). 

Just as a transfer function description of a system grows 
or shrinks with the number of inputs and outputs of the 
system, the structure matrix S also grows or shrinks with the 
number of system inputs and outputs. We call the number 
of outputs, p, the p''^-order resolution of the structural 
representation. Thus, when three of the states of a fourth 
order system are measured, the resolution of the structural 
representation is three, and Qt will be 3 x 3. The fourth state 
is hidden and does not appear in the third-order resolution of 
S in any way. Nevertheless, correct structural representations 
are consistent, in that zeros appear in lower-order (coarser) 
resolutions only if there are no hidden states from higher- 
order (finer) resolutions that could enable the interaction. 

Structure estimation for dynamic systems seeks to find 
S corresponding to a particular realization of a dynamic 
system, T, using only input output data. Before discussing 
how to solve this problem, however, we first outline two 
flawed approaches to this problem that appear from time- 
to-time in the literature. First, one may assume knowledge 
of the system order, n, and then proceed to attempt to infer 
information about structure in light of this assumption. Sec- 
ond, one may estimate a particular realization of T and then 
attempt to reconstruct S from the state space model. These 
approaches are not entirely unrelated, but we show next that 
either approach can easily lead to incorrect conclusions. 



A. Example: Erroneous System Order Assumption 

Although there are some reasonable techniques for esti- 
mating order from time-series data, there is no foolproof 
method available. In some applications, the most common 
technique for order estimation continues to be to assume 
that the measured outputs constitute the entire state vector, 
that is, that n — p. The following example demonstrates that 
making this assumption incorrectly may lead to completely 
erroneous structure estimates. 

Consider the network in Figure |l(a)| with three state 
variables structured in a chain, with the single input u driving 
xi, x\ in feedback with X3, and ^3 driving X2, characterized 
by the equations 
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From Xl and X2, we would like to be able to infer the structure 
u — > Xl — > X2, in spite of the fact that we may have 
no knowledge of the not-directly-observed, yet (indirectly) 
observable state X3. 
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(a) Network 



(b) Step response of X] and X2 



Fig. 1. Example of a simple 3 state network 

By assuming knowledge of the system order, one may then 
attempt to fit a state space realization directly from the data. 
In this case, any attempt to identify a second order system 
given the oscillating data shown above will result in an A 
matrix with complex eigenvalues. This implies that any real- 
valued A matrix that reasonably fits the data will have non- 
zero terms in its off-diagonal positions, leading incorrectly 
to a fully connected network structure estimate instead of the 
correct chain structure. ■ 

B. Example: Erroneous Structure from Realizations 

Suppose that after a sequence of experiments, one was 
able to identify the transfer function 
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It can be shown that this transfer function is consistent with 
two systems with very different structures, given by 
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The networks in Figure |2] correspond to each of the possible 
reaHzations of G. Note that without more information about 
the system, such that it is minimal, or order three, etc. then 
we would not be able to say anything about structure from 
the transfer function alone. 




(a) Possible network 
1. 



(b) Possible network 
2. 



Fig. 2. Two possible networks given the data. 

These examples demonstrate the difficulty of estimating 
network structure from data. Nevertheless, ideally one would 
estimate both the network structure and the system dynamics 
from data. In the next section, we introduce Dynamical 
Structure Functions as a mechanism for representing both 
system dynamics and structure. 

III. Dynamical Structure 

Consider the system given by Q. Given the special 
structure on C, we note that the first p state variables are 
actually the measured variables y. Defining x/, to be the 
remaining n — p "hidden" states, the system becomes 
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Taking Laplace Transforms of the signals, we then obtain 



Bi 
Bi 



U (8) 



sY ^ An An Y 
sXi, \ [ A21 A22 \ [ Xh 

From this equation it is easy to construct the transfer func- 
tions from the manifest variables z—[Y; U] to themselves. 
Solving for X/„ we have 



X,, = isI-A22r^A2iY+{sI-A22r'B2U 
Substituting into (jSjl then yields 

Y^WY + VU 

where W = An + A12 (i/ — A22)^' A21 and V = 
A12 A22)"' B2 + ^1- Let D be a matrix with the 
diagonal term of W, i.e. D = diag(Wi 1,1^22, ■■■,Wpp). Then, 

{si - D)Y ^ {W ~ D)Y + VU 



Note that W — D is a matrix with zeros on its diagonal. We 
then have 

Y = QY + PU (9) 



where 



and 



Q={sI-D)-^{W-D) 



(10) 

(11) 



The matrix g is a matrix of transfer functions from F, 
to Yj, i 7^ j, or relating each measured signal to all other 
measured signals (recall that Q is zero on the diagonal). The 
full transfer matrix from F to F thus becomes Qt — I + Q- 
Likewise, the transfer matrix from [/ to F is P. 

We thus can consider the transfer matrix, A^, relating all 
manifest variables, z, to themselves. This matrix is remi- 
niscent of the structure matrix S given in 0, except that 
the entries are transfer functions relating variables instead of 
binary values, this gives us the following definition. 

Definition 1: Given the system Q, we define the Dynam- 
ical Structure Function or Network, N, of the system to be 



N = 



I + Q P 
/ 



(12) 



where Q and P are as given in ilO\ and il 1\ . 

When this function is completely characterized by P and 
Q (when the system is open, that is, u represents completely 
free inputs unrelated to the measurements y), we refer to 
{P, Q) as the Dynamical Structure of the system. There are a 
number of properties of the Dynamical Structure Function 
that makes it useful for the structural analysis of linear 
systems: 

Proposition 1: Given the original realisation Q, every 
entry Nij is a strictly proper function and unique. 
Strict properness follows from the fact that (.?/ — D)^' (which 
is strictly proper) is multiplying transfer functions that are at 
most proper (never improper). This fact is important for the 
interpretation of as network structure. The directed edges 
associated with non-zero entries of this matrix imply causal 
relations; strict properness of the transfer functions preserve 
this interpretation. Uniqueness follows by construction of 
both Q and P. 

Proposition 2: The transfer function, G, of the system Q, 
is related to Dynamic Structure by 



G^iI-Q)-'P 



(13) 



This fact follows directly from ^ and F = GU and demon- 
strates that Dynamic Structure is a factorization of a transfer 
function into two parts, the output or internal structure, Q, 
and the input or control structure, P. 

It is now easy to see Nij = if and only if there is no direct 
or hidden connection from Zj to z,. The question is then on 
how to determine the — p and pm transfer functions in 
Q and P, respectively, to determine the Dynamical Struc- 
ture from data. This structure estimation, or reconstruction 
problem is addressed next. 



IV. Dynamical Structure Reconstruction 

Assume data is collected from the original system Q 
leading to the transfer function in Q relating Y = GU. 
Here we assume without loss of generality that G is full 
rank. Otherwise, there would be redundant inputs that could 
be removed to get a full rank G. Replacing Y = GU in 
equation and noting that the vector U is abitrarely yields 



by (|6} . In this case, equation has two equations and 
four unknowns 



{I-Q)G = P 



(14) 



This equation shows that there are more unknowns than 
equations and that in general Dynamical Structure of the 
p measurable states cannot be obtained from the m inputs. 
There are — p unknowns in Q, corresponding to all of the 
Qij which represents the internal Dynamical Structure. Then 
there are pm unknowns in P which represent the control 
Dynamical Structure on each measurable state. Thus, all 
together, there are a total of p^ — p + pm unknown but only 
a total of pm equations so the problem is under determined 
as we have p^ — p degrees of freedom. For instance, setting 
all Qij = (which means no connection between measured 
states) and f = G is a solution of but probably the wrong 
one. 

This clearly shows that the Dynamical Structure has more 
information than G and less than the original system Q. 
Thus, to find the Dynamical Structure from G we need more 
information. Either in the internal Dynamical Structure (if 
we know some Qij = 0), or on how the control is affecting 
measurable state (if some Pjj = 0). Next we assume we have 
no information on the internal Dynamical Structure (i.e. no 
information on Q) and consider the cases where: m < p (there 
are less inputs than measured states), m = p and m > p. 
Before that, we need the following technical result. 

Lemma 1: Rank(P) = rank(G). 

Proof: Since {I — Q)G = P, if suffices to show that 
rank(/- g) = p. It follows that rank(/- g) 



[l-{sI-D) ^ (An+An{sl-A22) 'Aji-d)} 



rank< 



= rank{s/-D- (Aii+Ai2(i/-A22) 'Aji^d)} 

= rank|i/- (Au +An{sI -AiiV^ Aji^} 

which has rank — p. ■ 

A. m < p: Less Inputs than Measured States 

lfm<p, i.e. there are less inputs than measured states, and 
we have no information on the internal Dynamical Structure 
then the Dynamical Structure cannot be recovered. To see 
this, note that in the best case scenario m = p — 1 and we 
would have mp = p^ — p equations. Since there are p^ — p 
unknowns from Q we would need to know P precisely. 

The example from section E-BI shows how different net- 
works satisfy il4\ if m < p. There we had two measurable 
states p = 2, a single input m=l and G = (Gii,G2i) given 
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We must solve for the internal Dynamical Structure {Qn 
and Q21) and the control Dynamical Structure {Pn and P2i)- 
Since there are only two equations, there are two degrees of 
freedom. A possible solution is to set Q12 = Q21 = 0, i.e. no 
internal connection between yi and y2- In that case. Pi 1 = Gn 
and P21 = G21 (Figure |2(a)t . Note that xt, is playing the role 
of a hidden state (as P21 is second-order) and the system 
is not controllable (due to the common pole at —1), which 
explains why G is second-order and there are three states. An 
alternative is to have P21 = (which fixes 221 = G2i/Gn) 
and Q12 = (which fixes Pn = Gn), which can be seen in 
Figure |2(b)| Note that in this case f 1 1 ^ as that would 
result in a non-proper Qi2- These two networks are different 
and obviously only one can be correct. 

B. m = p: Same number of Inputs than Measured States 

Theorem 1: If m = p and we have no information on the 
internal Dynamical Structure, then the Dynamical Structure 
can be reconstructed if and only if each input controls a 
measured state independently, i.e. P^j = for / ^ /'. In this 
case, the zeros of H — G^^ of the off diagonal define the 
internal Dynamical Structure and 

Hii 1 

nu tin 

Proof: The "if" part of the proof follows from the fact 
that there are p^ — p + p — p^ unknowns and p^ equations. A 
linear set of equations can be solved for Q and P. Multiplying 
on the right (/ - Q)G = P by H = G"' yields I Q = PH. 
Since Q has zeros on its diagonal and P is diagonal, we 
have 1 = PiiHii or Pa = 1 ///„■. Finally we can now solve for 
Q = I PH and the result follows. 

For the "only if", assume the Dynamical Structure can 
be reconstructed, i.e. we can solve for Q and P in il4\ 
uniquely and they are all strictly proper Since rank(G) = p, 
by lemma nrank(P) = p. Thus, there are at least p nonzero 
entries in P. 

To show that there are at the most p unknowns in P, 
assume there are additional unknowns in P, and there is 
a Dynamical Structure with strictly proper Q* and P* that 
satisfy J14t . We want to show that another Dynamical 
Structure with strictly proper Q ^ Q* and P ^ P* can 
be constructed. Consider a vector X stacked with all the 
unknown parameters, i.e. with all unknown Qij and Pij. 
Equation (I14> can then be written as s^X = SS, where both 
and ^ are functions of the elements of G. Because there 
are p'^ equations but 79^ + extra unknown elements in P, this 
system of equations is undetermined or jz/ has a null space. 
Let X 7^ be an element of the null space of ^ and X* 
contain the elements of Q* and P*, which satisfy .siX* — SS. 
Then, there exists a large enough positive integer n such that 

1 



is also a solution of s^X — and all the elements in X 
are strictly causal. We have then found another Dynamical 
Structure which contradicts the assumption. Thus, at the most 
there are only p unknowns in P. 

Finally, there must then be exactly p unknown and nonzero 
entries in P. Since P is full rank, each row and column must 
have exactly one of these entries. Without loss of generality 
the inputs can be renamed and reordered so that the diagonal 
of P contains the unknown and nonzero entries. ■ 

This theorem says that in addition to having a square and 
full rank G it is necessary and sufficient to know that each 
control / affects first state / before it affects any other measur- 
able state to reconstruct the Dynamical Structure. That allows 
to reduce the number of unknowns to — p + p = which 
can now be solved. 

However, if there is some a priori information about the 
internal Dynamical Structure (such as some of the Qij — 
0) then there is more flexibility and less information and 
constraints are required of Pij. As long as there are a total of 
p'^ nonzero elements between Qjj and Pjj then the Dynamical 
Structure can be reconstructed by solving the linear system 
of equations MA\ . 

If P is not diagonal and we have additional information 
on how the inputs affect the measured states, there may be a 
change of basis in the control vector that allows it to be 
converted to a diagonal matrix that can then be used in 
theorem n For example, if xi is controlled by ui+U2 and 
X2 by Ml — M2 then one could define two new input vectors 
v'l = Ml + M2 and v'2 = mi — M2. 

If all the states are measured and B = /, we have the 
following result. 

Corollary 1: If p = m = n and B — C = I then for / 7^ j, 
Hi] = Qij. Thus, Qij = (/ 7^ j) iff Hij = 0. 

Proof: The proof follows since is this case G{s) = 
(s/ — A)^' which means A. ■ 

Note that we if knew we were measuring all the states 
we did not need to know B. However, the information that 
the we measure all the state is not captured by ( I14t . unless 
we had imposed that there were only n modes available to 
construct Q and P. 

C. m> p: More Inputs than Measured States 

It may seem intuitive that if there are more inputs then 
there should be more information. However, the extra inputs 
are in fact redundant. The reason is the fact that although 
G is p X m, the rank(G) = p, which means that the inputs 
really only have m degrees of freedom. Thus, the problem 
reduces to having the same number of inputs as measured 
states. The difference here is that we may be able to choose 
from the m inputs p that are known to control directly each 
measurable state. 

D. The Danger of Steady-State Measurements 

Before ending this section, we want to clarify some 
misconceptions, especially some communities no so close 
mathematics, on steady-state identification versus time-series 
data. For instance, in [6] the authors proposed a method 



to estimate networks based on full state measurement and 
control. From corollary ^ the Dynamical Structure can be 
obtained from the zeros of the non-diagonal terms of H, 
which correspond directly to the entries A. However, in 
the realistic case there are less measurements and control 
available than states. If instead of estimating G from time- 
series data we were to use only steady-state data, this could 
lead to mistakes as the following example shows. 

Consider a third order system with measurements and 
control on the first 2 states xi and X2 and the following 
dynamics 

" -1 1 -1 
A= 1 -1 -1 
1 1 -1 

This is a fully connected network, and so we expect the re- 
duced network consisting on xi and X2 to be fully connected 
as well. In this case. 
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which could lead one to think the reduced order network 
is not connected at all, i.e. xi does not affect X2 and vice 
versa. In general, for third order systems this is always true 
if and only if 012033 +013032 ~ for the connection from 
X2 to xi and 021033 +023031 — for the connection from xi 
to X2- Note that even when these equalities are not exactly 
zero but near zero, the presence of noise may again lead to 
wrong decisions. 

V. Conclusions 

This paper discussed the role of network structure for LTI 
systems. In particular, it was shown that transfer functions 
alone contain no information about the structure of an LTI 
system. We then introduced a new representation for such 
systems, a factorization of the system's transfer function that 
we call Dynamical Structure. Dynamical Structure Functions 
contain more information about the system then the transfer 
function because they also describe the network structure 
between inputs and outputs. Nevertheless, Dynamical Struc- 
ture contains less information about the system than its state- 
space description because no attempt is made to realize the 
network structure relating the non-measured, hidden state 
variables to the rest of the system. In this way. Dynamical 
Structure is a convenient analysis tool somewhere in between 
a system's full state space realization and its transfer func- 
tion. 

We then used Dynamical Structure to explore the network 
reconstruction problem. In this problem, one would like to 
estimate network structure given only input-output data. This 
problem is extremely important for a variety of fields, such as 
biology or counter-terrorism, that attempt to draw structural 
conclusions from data. Necessary and sufficient conditions 
were presented that indicate that network reconstruction 



demands careful experiment design. Moreover, various ex- 
amples were provided throughout the paper that demonstrate 
how failure to respect the necessary conditions may lead to 
incorrect conclusions about the network structure. 
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